Di-jet hadron pair correlation in a hydrodynamical model with a quenching jet 
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In jet quenching, a hard QCD parton, before fragmenting into a jet of hadrons, deposits a frac- 
tion of its energy in the medium, leading to suppressed production of high-pT hadrons. Assuming 
that the deposited energy quickly thermalizes, we simulate the subsequent hydrodynamic evolution 
of the QGP fluid. Hydrodynamic evolution and subsequent particle emission depend on the jet 
trajectories. Azimuthal distribution of excess tt~ due to quenching jet, averaged over all the trajec- 
tories, reasonably well reproduce the di-hadron correlation as measured by the STAR and PHENIX 
collaboration in central and in peripheral Au+Au collisions. 
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From the general theoretical considerations, it was pre- 
dicted [l[ that in a dense deconfined medium, high-speed 
partons will suffer energy loss, significantly modifying the 
fragmentation function, which in turn will lead to sup- 
pressed production of hadrons. The phenomena called 
"jet quenching" , was later verified at Relativistic Heavy 
Ion Collider (RHIC), in Au+Au collisions at -^5^^=200 
GeV UMM- 

However it is unclear how the lost energy 
is transported in the dense medium. It has been sug- 
gested that a fraction of lost energy will go to collective 
excitation, call the "conical flow" p, H, 0|- The parton 
moves with speed of light, much greater than the speed 
of sound of the medium (c jet >> c s ), and the quench- 
ing jet can produce a shock wave with Mach cone an- 
gle, 6m = cos~ l c s /cj e t- Resulting conical flow will have 
characteristic peaks at 4> — tt ~ ®m and 4> = tt + 6m ■ 
Both in STAR % |f] and PHENIX [lfj experiments, in- 
dication of such peaks are seen in azimuthal distribution 
of secondaries associated with high px trigger in central 
Au+Au collisions. Mach like structure (splitting of away 
side peak) can also be obtained in various other models, 



e.g. gluon Cerenkov like radiation models [ll|,[12|], the 
parton cascade model fl3| , the Markovian parton scat- 
tering model Q, the color wake model Recently 
in [161 ] energy density wake produced by a heavy quark 
moving through a strongly coupled N=4 supersymmctric, 
Yang-Mills plasma is computed using ADS/CFT corre- 
spondence. Mach cone like structures is also observed 
for quark velocity greater than the speed of sound of the 
medium. 

Recently, we have numerically solved hydrodynamical 
equations with an (time dependent) source, representing 
the quenching jet 17|, [l8|, [l9( . It is assumed that the lost 
energy is quickly thermalised. Nevertheless, the energy is 
deposited locally along the trajectory of the fast parton, 
leading to local energy density inhomogencities, which 
if thermalised should in turn evolve hydrodynamically. 
Explicit simulation of hydrodynamic evolution with a 
quenching jet, indicate that the evolution of the fluid as 
well as subsequent particle emission are str ongly i nflu 
enced by the jet path length in the medium 
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Evolution produces 'distorted' Mach shock front like sur- 
faces. As predicted in [13], finite fluid velocity and also 
inhomogeneity of the fluid, distorts the Mach surfaces. 
Azimuthal distribution of pions, due to quenching jet 
only, do not show the characteristic of 'conical' flow. 
Rather, depending on the jet trajectory, it show a single 
peak either at cf> > tt or cf> < tt. The distribution averaged 
over all the jet trajectories, show two peaks with a dip at 
<t> = tt, mimicking the conical flow. However, in flij ]. we 
did not account for the limited px range (0.15 < px < 4) 
in the STAR experiment. Moreover, in [18(, jet trajec- 
tories were characterized by a single parameter, 4> pro d, 
while in a two-dimensional calculation, as in [l8| , unique 
characterization of jet trajectories need at least two pa- 
rameters. In the present paper, we have corrected the 
deficiencies in the model. Also, in addition to the STAR 
data H, we analyze the PHENIX data [10] on the di- 
hadron angular correlation in 0-5%, 5-10% and 60-90% 
centrality Au+Au collisions. 



A schematic representation of the jet moving through 
the medium is shown in Fig[TJ We assume that just be- 
fore hydrodynamics become applicable, a di-jet of high- 
Pt partons is produced. Strong jet quenching and sur- 
vival of the trigger jet, forbids production in the interior 
of the fireball. Jet pairs can be produced only on a thin 
shell on the surface of the fireball. For Au+Au colli- 
sions at impact parameter b, we assume that the di-jet is 
produced on the surface of the ellipsoid with minor and 
major axis, A = R - b/2 and B = R^/l - b 2 /4R 2 , with 
R = 6.4fm. One of the jet moves outward and escapes, 
forming the trigger jet. The other enters into the fire- 
ball. As shown in Fig[TJ the trajectory of the quenching 
jet can be uniquely characterized by two angles, 4> pro d, 
( tt < 4> P rod < tt) and (f> jet , {-tt/2 < 4> jet < tt/2). The 
fireball is expanding and cooling. The ingoing parton 
travels at the speed of light and loses energy in the fire- 
ball which thermalizes and acts as a source of energy 
and momentum for the fireball medium. We solve the 
energy-momentum conservation equation, 
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where the source is modeled as, 



(1) 



We solve EqsQ] in (r = \J1? — z^,x,y,ri 
coordinates, assuming boost-invariance. The source term 
Eqj3] violate the assumption of boost-invariance. Wc 
therefore modify it by replacing the (^-function in Q by 
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FIG. 1: Schematic representation of a jet moving through the 
medium. The high pr pair is assumed to produce at P on the 
surface of the fireball characterized by the angle 4> pro d- One 
of the jet escapes forming the trigger jet, the other move in 
the fireball at an angle <j)j e t- 

Massless partons have light-like 4-momentum, so the 
current J v describing the 4-momentum lost and de- 
posited in the medium by the fast parton is taken to be 
light-like, too. rjet(t) is the trajectory of the jet moving 
with speed \dx- }e t/dt\ =c. 4j| (x) is the energy loss rate 
of the parton as it moves through the liquid. It depends 
on the fluid's local rest frame particle density. Taking 
guidance from the phenomcnological analysis of parton 
energy loss observed in Au+Au collisions at RHIC [2l[ 
we take 
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where s(x) is the local entropy density without the jet. 
The measured suppression of high-p^ particle production 
in Au+Au collisions at RHIC was shown to be consis- 
tent with a parton energy loss of sLI| =14GeV/fm at 
a reference entropy density of so = 140 fm -3 [2l|. We do 
not considered the possibility that some jets might be 
stopped in the medium. It is implicitly assumed that 
the ingoing parton has enough energy to pass through 
the medium. The energy loss is weighted by the entropy 
density, as the fluid evolve, entropy density decreases and 
at the late stage of the evolution energy loss is minimum. 
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with a = 0.70 fm, rj_ = (x, y). 

Intuitively, this replaces the "needle" (jet) pushing 
through the medium at one point by a "knife" cutting 
the medium along its entire length along the beam direc- 
tion. Rather than 'conical', the replacement will produce 
a 'wedge' flow, over estimating the effect of jet quenching. 

The hydrodynamical equations are solved with the 
standard initialization described in 22], corresponding 
to a peak initial energy density of eo — 30 GeV/fm 3 at 
r = 0.6 fm/c. We use the equation of state EOS-Q de- 
scribed in [22j incorporating a first order phase transition 
and hadronic chemical freeze-out at a critical tempera- 
ture T c = 164 MeV. The hadronic sector of EOS-Q is soft 
with a squared speed of sound c 2 . ss 0.15. 

Some results of our simulations, for Au+Au collisions 
at impact parameter b=2.3 fm, are shown in FigJ5] In 
left panels (a),(b) and (c), for a few jet trajectories, con- 
tour plots of 'excess' energy density (energy density in 
evolution with a quenching jet minus the energy density 
in evolution without a quenching jet), after 8 fm of evo- 
lution, are shown. The contour plots of 'excess' energy 
density clearly show that the hydrodynamic evolution of 
the QGP fluid, in presence of a quenching jet, depend 
on the jet path length. For example, excess energy den- 
sity distribution in panel (a) and panel (c) are identical 
except for a rotation about w/2. For those two trajecto- 
ries, the jets traverse a similar path length in the fluid. 
We also note the excess energy density distribution in all 
the jet trajectories show Mach cone like surfaces. One 
also notices the distortion of Mach cone like surfaces. In 
panel (b) the excess energy density distribution is sym- 
metric with respect to the jet axis, but not in panel (a) 
or (c). Finite fluid velocity and inhomogen eity of the 
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medium distorts the Mach cone like surfaces 

In the right panels (d), (e) and (f) azimuthal distri- 
bution of 7r~ due to the quenching jet are shown. Us- 
ing the standard Cooper-Frey prescription, for each jet 
trajectory, we have calculated the azimuthal dependence 
of p T integrated (1 < p T < 2.5GeV) n~ yield (^r)jet 
at freeze-out temperature, Tp= 100 MeV. Azimuthal 
angle <fi is measured with respect to the quenching jet 
axis. We also calculate the n~ yield (^)nojet in an evo- 
lution with identical conditions but for the quenching 
jet. In (d)-(e) azimuthal distribution of excess pion yield 
(^rr)nojet, are shown. Excess ir~ distribution 
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do not show a two peak structure, rather it show a sin- 
gle peak either at > 7r or at < 7r. Only when the 
jet is along the diameter, one find a broad peak centered 
around tf> = iz. A similar picture is seen at other trajecto- 
ries also. Indeed, we find that, depending on the trajec- 
tory, a quenching jet produces peak either at 0_ ~ 2. ± 1 
rad or at </>+ ~ 4 + 1 rad. As it will be shown later, when 
many events are summed up, a two peak structure, akin 
to conical flow, appear. The conical flow or the splitting 
of the away side jet is an average effect, not to be seen 
in a single event. 
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FIG. 2: In left panels contour plot of excess energy den- 
sity due to a quenching jet in b=2.3 fm Au+Au collisions 
are shown. The di-jet is produced at <j> pro d = 0°, but with 
different orientations; (a) 4>j e ± = —45°, (b)(j>jet = 0° and 
{c)4>jet = 45° respectively. In the right panels, azimuthal 
distribution of tt~ , due the quenching jet only, from the cor- 
responding evolution are shown. 



In FigJ3J we have shown the perturbation produced in 
the momentum density AT 02 = T^ t (x,y) - T°l jet (x,y) 
due to a quenching jet. The trajectory of the jet is shown 
by the straight line. AT 02 do not remain confined along 
the jet trajectory, it deviates sideward. Even though the 
jet is restricted to move parallel to the x-axis and ex- 
pected to produce a peak at = 7r, due to finite fluid 
velocity and inhomogeneity, the momentum perturbation 
moves sideward and produce a peak at <p < n. We may 
mention that AT 02 rather look like that due to a de- 
flected jet In the deflected jet picture also , in a sin- 
gle jet event, the away jet is deflected sideward, resulting 
in a side ward peak. When many events are summed up, 
a double peak structure appears in the azimuthal distri- 
bution. Two-particle angular correlation can not discrim- 
inate between the two scenarios, e.g. the distorted Mach 



shock wave or the deflected jet. Three-particle angular 
correlation can possibly discriminate between them. 
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FIG. 3: contour plot of perturbation in momentum density 
T 02 (T 02 in evolution with a jet minus T 02 in evolution with- 
out a jet) in r — y plane, at a fixed x=2.6 fm. The trajectory 
of the jet is shown by the black line with arrow. Perturbation 
in T 02 do not remain confined along the jet trajectory. 



In FigHJ in 4 panels we have shown the di-hadron cor- 
relation as measured by the STAR || and the PHENIX 
[lfl ] collaboration. STAR collaboration measured the 
correlation function in 0-5% centrality Au+Au collisions. 
They measured charged hadrons ( 0.15 < pr < iGeV) 
associated with high pr (4 < p t ™" er < QGeV) trigger 
particle. The correlation function shown in Figd^, barely 
show the splitting of the away side peak. PHENIX collab- 
oration measured the correlation function as a function 
of centrality of collisions. They constructed the corre- 
lation function between the trigger (2.5 < pt < ^GeV) 
and charged hadrons (1 < p T < 2.5GeV). PHENIX cor- 
relation function in 0-5%, 5-10% and 60-90% centrality 
collisions are shown in Fig|4}),c and d. In 0-5% and 5- 
10% centrality collisions, the away side peak is splitted 
into two peaks, but not in 60-90% centrality collisions. 
In 60-90% centrality collisions show the usual structure 
in jet event, two peaks, one at <f> = and the other at 
<p = it. We also note that , splitting of the away side 
jet is more prominent in PHENIX than in STAR experi- 
ment, presumably due to more hard hadrons in PHENIX 
measurements. 

In FiglH solid lines are the pr integrated azimuthal 
distribution of excess pions in Au+Au collisions at im- 
pact parameter b=2.3 fm (panel (a) and (b)), b=4.1 fm 
(panel(c)) and b=12.1 fm (panel (d)) Au+Au collisions. 
They roughly corresponds to 0-5%, 5-10% and 60-80% 
centrality Au+Au collisions. The 7r~ production was av- 
eraged over all possible jet trajectories. For each jet tra- 
jectory {4> pro d, 4>jet), w e calculate the pt integrated ex- 
cess pions 't'jet) and average over all the pos- 
sible jet trajectories, 
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Azimuthal distribution of excess 7r~, normalized by a 
factor N norm w 2.5 — 3.5, reasonably well reproduces the 
shape of the di-hadron correlation in STAR and PHENIX 
measurements in central Au+Au collisions. It must be 
emphasized, the present model do not contain any pa- 
rameter, other than N norm . Considering that STAR and 
PHENIX measured charged hadron, of which only ~ 70% 
are charged pions, N norrn « 2.5-3.5, seems reasonable. 
The subtle difference between the STAR and PHENIX 
0-5% centrality data (splitting of the away side peak is 
prominent in PHENIX but not in STAR measurement) is 
also well reproduced in the model. The difference is due 
to pr range of associated particles. STAR data contain 
more soft hadrons than in PHENIX data. Hydrodynam- 
ics is not very successful in peripheral collisions. Still 
the model reproduces PHENIX data in 60-90% central- 
ity collisions. 
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FIG. 4: Di-hadron correlation as measured by the STAR @] 
and the PHENIX QU collaboration. The filled circles in (b)- 
(d) are (original) PHENIX data, the blank circles are gener- 
ated by reflection. The arrows in 4 panels indicate the approx- 
imate Mach peak position for average speed of sound 0.333 [7j. 
The solid lines are the present (hydro+jet) model calculation 
for the azimuthal distribution of excess n~ due to quench- 
ing jet, averaged over all possible quenching jets. They are 
normalized by a factor N norm . 

Apparently present results are at variance with the 
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conclusion of [17J. In [17j, it was concluded that the 
quenching jet possibly can explain the broadening of the 
away side peak but not the splitting. The conclusion 
was based on a single jet trajectory along the diameter. 
Distortion of Mach cone like surfaces in off-diameter tra- 



jectories produces the splitting. Mechanism of splitting 
of the away side peak in central and mid-central Au+Au 
collisions can now be understood in terms of jet quench- 
ing. An individual jet event does not produce 'conical' 
flow. Depending upon the trajectory, an individual jet 
produces associated particles either at <j)+ ± 8<f) > 7r or 
0_ ±6<fi < 7r. When averaged over all the jet trajectories, 
the associated particles show a two peaks structure, much 
akin to the 'conical' flow. However, we also notice that 
though the model correctly predict the splitting of the 
away side jet in to two peaks, the peaks in experiments 
are sharper than in the model calculations. The reason is 
possibly the assumption of boost-invariance. The 'needle' 
like jet is replaced by a 'knife', flattening the peaks. The 
model calculations also do not reproduce the exact the 
experimental peak positions. Simulated peaks are shifted 
from the experiment by about ~ 0.4 rad. As shown in 
[HI the peak in the azimuthal distribution of associated 
particles, depend, non-trivially, on the equation of state 
(speed of sound of the) medium. Deviation reflect our 
far from satisfactory knowledge of equation of state of 
the quarks matter for of the QGP/hadronic matter. 
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